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We examine the power spectrum of the energy level fluctuations of a family of critical power- 
law random banded matrices with properties similar to those of a disordered conductor at the 
' Anderson transition. It is shown both analytically and numerically that the Anderson transition is 

characterized by a power spectrum which presents l// 2 noise for small frequencies but 1// noise for 
■ larger frequencies. For weak diagonal disorder the analysis of the transition region between these 

two power-law limits provides with an accurate estimation of the Thouless energy of the system. As 
disorder increases the Thouless energy looses its meaning and the power spectrum presents a l// 2 
decay up to frequencies related to the Heisenberg time of the system. Finally we discuss under what 
circumstances these findings may be relevant in the context of non-random Hamiltonians. 
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I. INTRODUCTION 
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Level statistics is a powerful tool to investigate the properties of quantum complex systems. The spectrum is 
basis independent and typically easier to access either numerically or experimentally than the eigenfunctions. The 
analysis of the level statistics is usually carried out in two steps. First the spectrum is properly unfolded, namely, 
' by extracting the mean level density, the original spectrum is transformed into one with a mean level density equal 
to the unity. In a second stage the unfolded spectrum is analyzed by evaluating different spectral correlators. Two 
popular choices are the level spacing distribution P(s) (the probability of having two eigenvalues at a distance s) 
for short range correlations and the number variance S 2 (L) = (L 2 ) - (L) 2 (which measures the deviations of the 
number of eigenvalues in an interval L from its mean value) for long range correlations. In certain situations these 
t-H , spectral correlators present striking universal features. For instance, in the context of deterministic Hamiltonians the 
■ celebrated Bohigas-Giannoni-Schmit conjecture [| states the level statistics of a deterministic quantum system whose 
| classical counterpart is fully chaotic depends not on the microscopic details of the Hamiltonians but only on the global 
symmetries of the system and are identical to those of a random matrix with the same symmetry, usually referred to 
as Wigner-Dyson statistics (WD) . 

Remarkably the same WD statistics also describes y| the spectral correlations of a disordered system in the 
' metallic limit. In the strong disorder limit localization sets in, the spectrum is not correlated and the level statistics 
is universally described by Poisson statistics. For deterministic systems the same statistics is generic of systems 
whose classical dynamics is integrable £| . From a practical point of view, universality is usually tested by comparing 
correlators like P(s) or X 2 (£) of a specific system with the universal predictions of WD or Poisson statistics. 

Recently a different spectral characterization based on techniques of time series was introduced in the context of 
quantum chaos 0,0- In 0,0 the unfolded energy spectrum is formally considered as a discrete signal and the energy 
levels are interpreted as a time series. Specifically they investigate the power spectrum, 
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S n = ^(si- s) = e n+ i - ei - n, (2) 
i=i 

where Si = ej+i — e*, is ith unfolded eigenvalue and N is the length of the series. The correlator S n thus gives the 
deviation of the ith nearest neighbor spacing Si from its mean value s which by definition is the unity for unfolded 
eigenvalues. 

It was found that the S(k) for systems whose classical counterpart is completely chaotic or integrable has a universal 
power-law form with an exponent depending on the classical dynamics: S(k) ~ 1/k for chaotic and S(k) ~ 1/fc 2 for 
integrable motion. 
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Universality in the spectral correlations has also a counterpart in the eigenfunctions properties. Thus Poisson 
statistics is associated with exponential localization of the eigenfunctions and WD statistics is typical of systems in 
which the eigenstates are delocalized through the sample and can be effectively represented by a superposition of 
plane waves with random phases. 

Despite its robustness, these universal features are restricted to long time scales (related to energy scales of the 
order of the mean level spacing) such that an initially localized wave-packet has already explored the whole phase 
space available. In other words, universality is related to certain ergodic limit of the quantum dynamics. For shorter 
time scales the system has not yet relaxed to the ergodic limit and deviations from universality are expected. For finite 
disordered systems this scale is given by the dimensionless conductance g = E c j A (E c , the Thouless energy, is a scale 
of energy associated with the classical diffusion time through sample and A is the mean level spacing) which roughly 
speaking gives the number of eigenvalues whose spectral correlations are universally described by WD statistics. 

Striking universal features not related to any ergodic limit (they persist beyond the mean level spacing scale) have 
also been observed in a disordered system at the metal-insulator transition also referred to as Anderson transition 
(AT). It is by now well established that a disordered system with short range hopping in more than two dimensions 
undergoes an AT 0, Q at the center of the band for a critical amount of disorder (for critical we mean a disorder 
such that if increased all the states in the band become exponentially localized) . Since the dimensionless conductance 
g ~ 1 is about the unity at the 3D (or 4D) AT, the level statistics for eigenvalues separations larger than the mean 
level spacing describes truly dynamical features of the system. By universality we mean that the level statistics do 
not depend on boundary conditions, shape of the system or the microscopic details of the disordered potential though 
some features may depend on the dimensionality of the space. 

Signatures of the AT are found in both the level statistics and the eigenfunctions. Systems belonging to this new 
universality class have multifractal eigenstates. Intuitively multifractality means that the eigenstates have structures 
at all scales. In a more formal way multifractality is defined through the anomalous scaling of the eigenfunctions 
moments V q = J d d r\tp(r)\ 2 ' 1 with respect to the sample size L as V q oc L-^fa-i), where D q is a set of different 
exponents describing the AT Q. 

Level statistics at the AT (commonly referred to as 'critical statistics 10]) is intermediate between WD and Poisson 
statistics. Although a formal definition is still missing, typical features of critical statistics include: scale invariant 
spectrum level repulsion and linear number variance (^(L) ~ \L) j3 as for a insulator (x = 1) but with a slope 
X < 1 (0.27 for the 3D Anderson transition). Similar spectral properties has also been found in random matrix models 
based on soft confining potentials ^| , effective eigenvalue distributions [Til PHI related to the Calogero-Sutherland 
model 

M at finite temperature and random banded matrices with power-law decay 0. The latter is specially 
interesting since an AT (for the case 1 fr decay) has been analytically established by mapping the problem onto a non 
linear a model. 

In this paper we propose an alternative spectral characterization of the Anderson transition based on the analysis of 
the power spectrum S(k) introduced above. We shall also see this technique provides with an accurate way to locate 
the Thouless energy of a disordered system. Finally we will discuss the relevance of our findings in the context of non 
random system. It will shown that S(k) ~ 1/k 2 is not always related to integrable classical motion. Consequently a 
precise classificatory scheme based on S(k) must have into account other features of S(k) besides the exponent of the 
power-law decay. 

The organization of the paper is as follows. In the next section the model to be investigated is introduced. The 
power spectrum S(k) is evaluated both analytically and numerically for a broad range of parameters in section three 
and four. From these results we present a novel spectral characterization of the Anderson transition based on the 
analysis of S(k). We will also argue that this technique can be utilized to detect the Thouless energy in a disordered 
system. Finally in section five we discuss in what situations our findings may be relevant for non-random quantum 
systems. 

II. THE MODEL 

In this section we evaluated both analytically and numerically the power spectrum S(k) of a critical random banded 
matrix associated to critical statistics. 

Unlike WD or Poisson statistics, critical statistics is not parameter free. Together with universal features such as 
scale invariance, level repulsion and linear number variance S 2 (L) = xL i » 1, there are also system dependent 
features as the numerical value of the slope of the number variance x- For instance, for short range Anderson models, 
this slope depends on the Euclidean dimension d of the sample. Thus for the lower critical dimension d = 2 + e (e > 0) 
§,X~e<l- In the opposite limit d 3> 1, x < 1 is close to the unity similar to the case of an insulator. 

In this letter, instead of studying directly these short range Anderson model, we will focus on certain generalized 
random matrix models which has been shown to reproduce critical statistics [l(| with great accuracy. An advantage 
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FIG. 1: Power spectrum S(k) as a function of k. Symbols represent the numerical results (the matrix size is 3000 and the 
number of eigenvalues considered is N = 1024 around the center of the band) for the critical random banded model Eq[3] 
for different bandwidth b, the power spectrum was evaluated from Eq. Lines represent the analytical prediction of critical 
statistics Eq. Illl with the TLCF given by Eq. |S| For all b the agreement between theory and numerics is excellent. As predicted 
for 6 > 1 we observe two different regions: N/k ^> 2ixb corresponding with S(k) ~ l/k' 2 and N/k <C 2nb with S(k) ~ l/k 
similar to the prediction of WD (GUE). However for b <C 1, S(k) ~ l/k 2 for almost all k. 



of these models is that exact analytical solutions are available in a certain region of parameters 0, 0, . 

We investigate the ensemble of random complex Hermitian matrices H . The matrix elements Hij are independently 



distributed Gaussian variables with zero mean (Hi, 



and variance 



1 sin 2 (7r(i- j)/N) 
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For any value of the bandwidth < b < oo, the spectral correlations are given by critical statistics and the eigenvectors 
are multifractal exactly as at the conventional Anderson transition in 2 < d < oo |17| . The limit b — * oo corresponds 
with the standard Gaussian Unitary Ensemble (GUE) of random matrices. The region b 3> 1 (weak diagonal disorder, 
X <C 1) corresponds with d = 2 + e (e <C 1) and the 6 <§; 1 limit with d 3> 1 and y < 1 (strong diagonal disorder). For 
Hermitian matrices these two limits are accessible to analytical techniques^! ha. Here we do not discuss the details 
of these calculations but just enumerate certain results we will use later on in the calculation of the power spectrum 
S(k). 

For i» > 1 the level statistics can be rigorously investigated after mapping the random banded matrix onto a 
supersymmetry sigma model. It can be shown 17] that in this limit the connected part of two level correlation 
function (TLCF) is given by, 



(4) 
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where p(s = E/A) = J^i $( s i ~ s ) i s the spectral density in units of the mean level spacing A = l/(p(0)) and brackets 
stand for ensemble average. 

For b <C 1, various spectral correlators can also by calculated explicitly by a recently developed virial expansion 
around the Poisson limit (f° r a more heuristic approach see 01 )• 

The intermediate re gion of b ~ 1 is not yet accessible to analytical techniques. However there is a closely related 
random matrix model pj| which is exactly solvable for any b and which has the same TLCF (to leading order in b) 
in the two regions (b <C 1. b 3> 1) above discussed. Its joint probability distribution is given by 

P(H,b) = J dUe -h™ ] e -^V,HWM\ (5) 

Here, the N x N matrices H and U are Hcrmitian and Unitary, respectively, and the integration measure dU is the 
Haar measure. Despite its complicated form, it can be shown [141 [l5j that the joint distribution of eigenvalues of 
H is equal to the diagonal element of the density matrix of a system of free spinless fermions at finite temperature 
T = l/b confined in an harmonic potential. By using elementary statistical mechanics techniques it can be shown 
that the TLCF for arbitrary b is given by, 

cos(7rsi/p(0)) 1 
7(0) lTze* 



Ra(») = 6(a) - ( / ""^^ r—$ dt ) (6) 



o 



where h = \/2i:b,z = e i/l_ 1 and p(0) = 1+ l z , t -i ■ We shall use this expression in the analytical evaluation of the 
power spectrum for intermediate values of b and then check its validity by carrying out numerical simulations of the 
random banded model Eq. 




FIG. 2: Power spectrum S(k) versus N/k in the limit b < 1. Crosses and circles represent the numerical results (the matrix 
size is 3000 and the number of eigenvalues considered is N = 1024 around the center of the band) for the critical random 
banded model Eq. I^Jfor b = 0.25, 1 respectively. The power spectrum was evaluated from Eq. |7| The dashed line corresponds 
to the best fit S(k) — l/k a (the error in a is about Ao = ±0.02) in the limit N/k 3> 2nb and the solid line is the prediction 
of WD statistics obtained from the GUE. The intersection between the linear fit and the WD statistics prediction corresponds 
in principle with the Thouless energy of the system. However we observe that, for b < 1, S(k) is always different from the WD 
prediction and consequently no clear Thouless energy can be defined. 



III. ANALYTICAL EVALUATION OF S(k). POWER SPECTRUM CHARACTERIZATION OF 

CRITICAL STATISTICS. 



Our goal is to compute the power spectrum 

(S(k)) = ( 



1 N 



exp 



-2nikn 



(7) 
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with S n = X)"=i s i~ s = e n+i — £i — n, e n represents the nth unfolded eigenvalue of the critical random banded model 
Eq. © and the brackets stand for ensemble average. 

In a first stage we evaluate (S(k)) in a continuous approximation, (S(k)) — (S(k)) cont , ~* J \ s i ^ p( £ )de and 
s — > p(e')de' where p(e) — ^ 5(e — e^) is the full spectral density and p(e) is just the smooth part of it (the one 
utilized to unfold the spectrum). 

The power spectrum is now given by, 



(S(k)) c 



deS f (e) exp (— 2-Kek) 



' /-+oo r-\-oo /*e' /•€ \ 

/ / de'decxp(-2iri{e- e')k) / / p(a)p(a')da'da) (8) 

1 J — OO J — OO J J I 



with J'(e) = J e dep(e) and /?(e) = p(e) — p(e). 

After integrating by parts in e and e' the above expression simplifies to, 
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(S(fe)) cont = ^2 y y de'deexp(-27r l (e-e')fc)i? 2 (e,e') (9) 

where i?2( e ! e') is the TLCF defined previously. Since the spectrum is translational invariant (if we are far from the 
edges) i?2(e, e') = #2(3 = e — e') and 

1 r 00 1 
= i^2p /_ exp(-27rz S fc)i? 2 ( S ) = -^K(k) (10) 

where K(k), the Fourier transform of the TLCF, is usually referred to as the spectral form factor. 

Once we have obtained an explicit expression for the power spectrum in terms of known quantities as i? 2 (s) we 
have to go back to the original discrete formulation. This can be easily done by following standard relations between 
the discrete and the continuous Fourier transform, we only present the final result and refer to 0, l20| for additional 
details, 

K $) ■ >fi K{t + q) , K(q-t) , Q 

with t = k/N and (3 a constant given by (3 — [K(0) — 1] 2 /12 which accounts for the differences between the fluctuations 

of J e f e (p(a)p(a'))da'da and those of the original discrete correlator S n (see |2(| for details). For the sake of simplicity 
we set (S(k)) — S(k). The above expression combined with Eg 14161 provides with a closed and compact expression for 
the power spectrum as a function of the band size b = l/2nh. 

In the region of b 3> 1 the spectral form factor can be explicitly evaluated by using the TLCF of Eq. . 



K(t) = l-\ 



2 — 2 A /2 + 2i 

(1 - t) coth ( — - — + (1 + 1) coth — - — I - It coth {It jh) 



(12) 



We can distinguish two different regions. For t <C h = 1/2tt&, corresponding to eigenvalues separated a distance 
much larger than the mean level spacing, K(t) ~ h/2 is a constant and S(k) ~ h/(8ir 2 t 2 ), similar to the case of 
Poisson statistics. However for Poisson K(0) = 1 but in our case K(Q) = h/2. This is an important difference since 
Poisson statistics is associated with eigenstates exponentially localized but for K(0) — \ 7^ 1 (where x 1S the slope of 
the number variance) the eigenstates are multifractal. 

It seems that, at least in this case, the exponent of the decay of S(k) does not completely specify the nature of the 
quantum motion. We will go back to this point when we discuss applications of our work in the context of quantum 
chaos. 

In the opposite limit t ^> h — 1/2-Kb, Kit) = t and S(t) ~ l/(27r 2 £) in agreement with the result for WD statistics 
(GUE). 

The transition region separating the two types of decay (1/f and 1/t 2 ) corresponds to the Thouless energy of the 
system. As usual it separates short range correlations still controlled by WD statistics from larger scales in which 
typical features of the AT appear. 

We have thus found that the power spectrum in the limit b 3> 1 corresponding to the case of a disordered system in 
2 + e dimensions with short range disorder at the AT has different power-law decay depending of the spectral region 
of interest, these differences can be effectively utilized to find signatures of an AT from a given spectrum. 

Analogously, in the region b = l/2nh <C 1, which corresponds with the case of disordered conductor in d ^> 1, a 
straightforward calculation shows that K(t) = 1 — 1/ y2h for t <C h and then goes to K(t) = 1 for t ^> h. Consequently 
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S(t) ~ 1/t 2 up to scales smaller than the mean level spacing. Strictly speaking there is a narrow transition region 
already for t > 1 in which K(t) is linear However it is difficult to interpret it as a Thoulcss energy since even for scales 
shorter than the mean level spacing the spectral correlations are different from WD statistics. 

Finally we mention that for intermediate b there is not a rigorous analytical relation for the TLCF. However we 
shall see that the conjecture Eq © describes very accurately the numerical results (see Fig 1). 




FIG. 3: Power spectrum S(k) versus N/k in the limit b 3> 1. Crosses and circles represent the numerical results (the matrix 
size is 3000 and the number of eigenvalues considered is N = 1024 around the center of the band) for the critical random 
banded model Eq[3] for 6 = 4,16 respectively. The power spectrum was evaluated from Eq. |7| The dashed line corresponds 
to the best fit S(k) = l/k a (the error in a is about Aa — ±0.02) in the limit N/k 3> 2nb and the solid line is the prediction 
of WD statistics obtained from the GUE. The intersection between the linear fit and the WD prediction corresponds with the 
Thouless energy of the system. The units has been chosen such that the intersection point gives the dimensionless conductance 
of the system g. Thus g ~ 25 for b — 4 and g ~ 100 for b — 16. 



IV. NUMERICAL CALCULATIONS 



We now investigate numerically the random banded model in Eq|3]in order to test the analytical predictions of the 
previous section. 

We studied the power spectrum S(k) by direct diagonalization of the critical random banded matrix Eq[3] for 
different matrix sizes (almost all of our plots are for 3000 though we try higher volumes, up to 5000, in order to check 
that our results are not size dependent). The number of different realizations of disorder is chosen such that for each 
matrix size the total number of eigenvalues be at least 2 x 10 6 . Typically around 35% of the eigenvalues around the 
center of the band are utilized. The eigenvalues thus obtained are unfolded with respect to the mean spectral density. 
The power spectrum is calculated by using Eq0 where the Fourier transform is evaluated by using a fast Fourier 
transformation routine. 

In Fig 1 we have plotted log S(k) for different b. In all cases the matrix size was 3000 and the evaluation of S(k) was 
carried out within a band around the center of the spectrum containing 1024 eigenvalues. As observed the agreement 
between the analytical (with J?a(s) given by Eq.© and h = l/2nb) and numerical results is excellent for all b. Also 
in agreement with the analytical prediction we observe that, for 6> 1, the power spectrum switches from S(k) ~ 1/k 
for N/K <C h — 1/2-7T& to S(k) ~ 1/k 2 in the opposite limit. However for k< 1, S(k) ~ 1/k 2 for almost all accessible 
frequencies. 

We conclude after the analytical and numerical analysis that the AT in a disordered conductor can be satisfactorily 
detected and examined by looking at the power spectrum of a signal 5 n consisting of the fluctuations around its mean 
value of the nearest neighboring spacings S{ = ef+j — £j. 

We have also found that S(k) provides with an accurate method to locate the Thouless energy of a generic disordered 
conductor. As mentioned previously, the Thouless energy E c is a scale of energy related with the classical diffusion 
time through the sample. In units of the mean level spacing A it gives the dimensionless conductance g, namely, 
the number of eigenvalues g = E c / A for which the universal results of WD statistics apply. From a practical point 
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of view the evaluation of g from a given spectrum is a hard task since it may depend on what spectral correlator is 
used. Thus S 2 (L) gives a prediction of g bigger than that of P(s) but much smaller than that of the spectral rigidity 
A3 (L) (see for a definition) . Another problem is that even for each particular correlator the value of g is somewhat 
ambiguous since it is far from clear how to locate even approximately the point in which WD ceases to be applicable. 
Below we show that S(k) provides with a more efficient and precise way to locate and analyze g. The idea (see Fig 2 
and Fig 3) is to plot S(k) as a function of N/k. Then in the region N/k ^> 2irb we fit S(k) to a linear (in a log scale) 
curve S(k) ~ l/k a (from the previous analysis a ~ 2). In the opposite limit S(k) should be given by the prediction 
of WD statistics. We define the Thouless energy as the intersection between the WD prediction and the linear fit. 

In Fig 2 we see that for b ^> 1 the intersection of these two curves gives g ~ 2nb 3> 1 in good agreement with the 
theoretical prediction. However in the region b < 1 (see Fig 3), though formally a Thouless energy can be defined 
thorough the intersection of the two curves, its interpretation as the limit of applicability of WD statistics is dubious 
since even in the limit N/k -C 2irb deviations with respect to the WD prediction are clearly visible. 

Finally we mention that, as observed in Fig 1, the best fit of the numerical results does not occur at the analytical 
estimation h = l/2nb. There are two reasons for that disagreement: The analytical results are strictly valid only at 
the center of the band. Eigenstates beyond this region are still critical but are described by an effective bandwidth 
|19|| smaller than b. On the other hand finite size effects are important in the limit t — > since we are testing the 
largest eigenvalue separations. However we have decided not to reduce the spectral window in order to give a full 
global picture of the power spectrum at the AT. After all, as shown in Fig 1, these effects are easily compensated by 
slightly modifying b. 




FIG. 4: (Left) Power spectrum S(k) obtained from 25 set of 256 eigenvalues after the numerical diagonalization of the evolution 
matrix associated to Eq. 11311 with V(q) = 0.2Ln\q\ (cross). The dashed line is the prediction of SP statistics, the dotted-dashed 
line is the best fit S(k) ~ l/k a and the dotted line is the prediction of WD statistics as obtained from the Gaussian Orthogonal 
Ensemble (GOE) of random matrices. (Right) Poincare section from a single initial condition po = 0.2 and go = 0.6 after 
300000 iterations. 



V. APPLICATION TO QUANTUM CHAOS 



In this final section we investigate possible applications of our previous findings in the context of quantum chaos 
and also argue that S(k) ~ 1/k 2 is not exclusive of quantum deterministic system whose classical counterpart is 
integrable. 

Critical statistics and multifractal wavefunctions are very much universal so they should also appear in deterministic 
quantum systems. Indeed, in a recent letter [2lT | we have established a novel relation between the presence of anomalous 
diffusion in the classical dynamics, the singularities of a classically chaotic potential and the power-law localization 
of the quantum eigenstates. Specifically, it was found that for a kicked rotor with Hamiltonian 



n = p - + v{ q )Y J 5{t-nT) 

n 



(13) 
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(with q e [—1, 1)) both level statistics and eigenfunctions are similar to the ones at the AT (critical statistics) pro vided 
that V(q) has a log (in the simplest case V(q) = e Ln\x\) or a step-like singularity. It was also found in |2l| that 
these findings are universal in the sense that neither the classical nor the quantum properties depend on the details 
of the potential but only on the type of singularity. Deviations from WD statistics not coming from a mixed phase 
space has also been reported in a variety of systems: Coulomb billiard |2^|, Anisotropic Kepler problem [24^. a kicked 
rotor in a well potential [2£j and pseudointegrable billiards [2H Hl\ . For the latter it was found (2|| that the level 
statistics is accurately described by a the classical Dyson gas with the logarithmic pairwise interaction restricted to 
a finite number k of nearest neighbors. Analytical solutions are available for general k. For k = 2, usually referred 
to as semi Poisson (SP) statistics, R 2 (s) = 1 - e~ 4s , P(s) = 4se~ 2s and E 2 (L) = L/2 + (1 - e~ 4i )/8. We have also 
found [2^ that SP statistics describes accurately the spectral correlations of the above kicked rotor with a step like 
singularity and also provides with a reasonable description for the log singularity but only for e ~ 0.2. 
Due to the simplicity of the TLCF in SP statistics one can evaluate the power spectrum exactly, 



S(k) ' 



1 



16 + AttH 2 



(14) 



Thus the power spectrum associated to SP statistics has 1/t 2 decay even though the classical dynamics is not 
integrable, this is also in agreement with the prediction of critical statistics for b <C 1. Thus S(k) ~ 1/t 2 (t = k/N) 
is not always a signature of classical integrable dynamics. Although generically the power spectrum associated with 
classically integrable systems has this feature, other types of non integrable dynamics may have 1/t 2 as well. In order 
to fully characterize the classical dynamics from S(k) one has to specify not only the exponent but also and additional 
point of the curve, for instance K(0). The point is that a 1/t 2 decay only tell us that the form factor is constant. 
However, as mentioned previously, the physical properties of the system are strongly modified by a spectral form 
factor different from the unity. 

As a further corroboration of our claims we have evaluated numerically S(k) for the Hamiltonian of Eq. 1131 with a 
potential V(x) — Q.2Ln\x\. We diagonalize numerically the evolution matrix associated to the Hamiltonian Eg 1 131 for 
N = 5200, S(k) is obtained from Eq0 In order to improve statistics we divide the original spectrum in 20 set of 256 
eigenvalues. As shown in Fig 4, S(k) ~ 1/t 2 for almost all t in close agreement with the prediction of SP or critical 
statistics. However (see right plot) the associated Poincare section obtained from just a single initial condition is very 
different from that of an system whose classical dynamics is integrable. 



VI. CONCLUSIONS 



We have shown that the power spectrum of the energy level fluctuations at the Anderson transition for a family of 
critical power-law random banded matrices is characterized by a power spectrum which l// 2 noise for small frequencies 
and 1// noise for larger frequencies. In the weak disorder limit b 3> 1, the analysis of the transition region between 
these two power-law limits provides with an accurate estimation of the Thouless energy of the system. As disorder 
increases the Thouless energy looses its meaning and the power spectrum presents a l// 2 noise up to frequencies 
related to the Heisenberg time of the system. Finally we discuss under what circumstances these findings may be 
relevant in the context of non-random Hamiltonians. Specifically it is shown that the exponent of the power-law decay 
of S(k) does not fully specify the type of motion of the classical counterpart. 
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